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Abstract 



Let s, t, m, n be positive integers such that sm = tn. Let M(m, s; n, t) be the 
number ofmxn matrices over {0, 1,2,...} with each row summing to s and each 
column summing to t. Equivalently, M(m, s;n,t) counts 2- way contingency tables 
of order m x n such that the row marginal sums are all s and the column marginal 
sums are all t. A third equivalent description is that M(m, s; n, t) is the number 
of semiregular labelled bipartite multigraphs with m vertices of degree s and n 
vertices of degree t. When m = n and s = t such matrices are also referred to as 
n x n magic or semimagic squares with line sums equal to t. We prove a precise 
asymptotic formula for M(m, s; n, t) which is valid over a range of (m, s; n, t) in 
which m, n — > oo while remaining approximately equal and the average entry is not 
too small. This range includes the case where m/n, n/m, s/n and t/m are bounded 
from below. 

1 Introduction 

Let m, s, n, t be positive integers such that ms = nt. Let M(m, s; n, t) be the number of 
m x n matrices over {0, 1,2,...} with each row summing to s and each column summing 
to t. Figure [1] shows an example. 
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5 4 11 
2 8 9 1 
8 3 6 3 



20 
20 
20 



15 15 15 15 



Figure 1: Example of a table counted by M(3, 20; 4, 



15). 



The matrices counted by M(m, s; n, t) arise frequently in many areas of mathematics, 
for example enumeration of permutations with respect to descents, symmetric function 
theory, and statistics. The last field in particular has an extensive literature in which such 
matrices are studied as contingency tables or frequency tables; see [17] as an example of 
enumerative work, [22] for a lengthy bibliography, and [H] for a survey of applications. 
The matrices counted by M(m, s;n,t) have also regularly been the topic of papers whose 
primary focus is algorithmic, both with respect to the problem of generating contingency 
tables with prescribed margins at random, and with respect to approximate counting. For 
random generation see, for example [13 [31]; for approximate counting see the recent [5]. 
Other important recent studies of both random sampling and approximate counting in- 
clude HHEI]. 

A historically early study of integer matrices with specified row and column sums 
appeared in MacMahon [25J, where the numbers M(m, s; n, t) are identified as coefficients 
when certain functions are expanded in terms of standard bases of symmetric functions. 
In a later paper [26], MacMahon studied "general magic squares" — n x n integer matrices 
with all row and column sums equal to a prescribed value t. Stanley [3H Chap. 4] refers 
to these as magic squares with line sums equal to t, and Stanley's [321 Chap. 1] provides 
further history of this topic. For fixed n, M(n, t; n, t) is a polynomial in t known as the 
Ehrhart polynomial of the Birkhoff polytope. This has been computed in closed form for 
n < 9 [2J. (The Birkhoff polytope is the set of all doubly stochastic real matrices.) The 
results in the present paper enable us to give the first asymptotic formula for the volume 
of this famous polytope [9]. 

Our focus in this paper is the asymptotic value of M(m, s; n,t). The first significant 
result on asymptotics was that of Read [33], who obtained the asymptotic behavior for 
s = t = 3. In the four year period from 1971 to 1974 three published papers [SHU [16] gave 
the asymptotic value for the case s, t bounded. To our knowledge, this was the extent 
of published, rigorously proven formulas until the present work. In a paper currently 
in preparation [2TJ, Greenhill and McKay obtain a formula valid in the sparse range 
st = o((mn) 1 ' 2 ) . 

Define the density A = s/n = t/m to be the average entry in the matrix. As early 
as [HI p. 100], Good implicitly gave the estimate M(m, s; n, t) ps G(m, s; n, t), where 
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G(m, s; n, t) 



n + s — 1\ fm + t 
s J [ t 



mn + Xmn — 
Xmn 

Good spelt this out explicitly in [19J and again with Crook in [20] ■ In the first paper, Good 
gave a heuristic argument for this approximation based on steepest-descent, "leaving aside 
finer points of rigor." His calculation also yielded another approximation [19, Eqn. B2.22] 
that is very similar when m « n but larger otherwise. We will show in this paper 
that G(m, s;n,t) is a remarkably accurate approximation of M(m,s;n,t), being out by 
a constant factor over a wide range and perhaps always. Another important estimate of 
M(m, s;n,t) was developed by Diaconis and Efron [15] . It aims for accuracy when the 
sum of the matrix entries Xmn is very large, and is quite poor when that condition is not 
met. 

Our main result is as follows. 

Theorem 1. Let s = s(m,n),t = t(m,n) be positive integers satisfying ms = nt. Define 
X = s/n = t/m. Let a, b > be constants such that a + b < |. Suppose that m, n — > oo in 
such a way that 

' I + — + —)< alogn. (1.1) 



4A(1 + A)V Qn Qm 
Then 

M(m,s;n,t) = G(m, s; n, t) exp(± + 0(n' b )) (1.2) 

(A- A (l + A) 1+A ) mw 

~ (47rA) (m+n_1)/2 m (n " 1)/2 n (m " 1)/2 

(\ l + 2A/m n\ h ,\ 



Corollary 1. Under the conditions of Theorem^ if in addition mn/X 2 — > ; 

M(m, S ;n,t) = (A + |) { ^ 1)( "- 1) ^^exp(I + 0(mn/A 2 + n- 6 )). 



m\ n\ 

Remark 1. The following is inspired by a similar observation made by Good /TP) /. The 
number of k-tuples of nonnegative integers satisfying n x + • ■ ■ + n k = K is well 
known to be (^^j^ 1 ) = ( K+ k~ 1 ) • This allows an instructive interpretation of ( 11. 2ft as 
M(m, s; n, t) = NP^E, where 

N= /mn + Xmn-l\ _ x /„ + a - A™ ^(m + t 

\ Xmn J \ s J \ t 

E = E(m, s; n, t) = exp (| + 0(n~ b )) . 
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Clearly N is the number of tables whose entries sum to Xmn = sm = tn. In the uniform 
probability space on these N tables, P 1 is the probability of the event that all the row sums 
are equal to s, and P 2 is the probability of the event that all the column sums are equal 
to t. The final quantity E is thus a correction to account for the non-independence of 
these two events. 

In an earlier paper [S] we showed how B(m, s; n, t), the number of such matrices whose 
entries are taken from {0, 1}, can be computed exactly for small m,n. That algorithm is 
easily adaptable for efficient computation of exact values of M(m, s; n, t), and the values 
so computed suggest the following conjecture. 

Conjecture 1. Consider a 4-tuple of positive integers m, s, n, t such that ms = nt. Define 
A(m, s; n, t) by 



,'m + 1 \ {ra-l)/2 fn _|_ 1 ^ („-l)/2 

M(m, s;n,t) = &(m, s;n,t) 



m / \ n 



( i A(m, s:n, t) 
xexp -1 + 

\ z m + n 

Then < A(m,s;n,t) < 2. 



;i-4) 



The factor (1 + l/m)*-" 1-1 ^ 2 , which approaches e 1 ^ 2 as m — > oo, and the similar factor 
(1 + l/ra) 1 -™ -1 ^ 2 , appear naturally in the analysis of M(m, s;n,t) when one of m,n goes 
to oo much faster than the other. This can be seen in [17, eqn. (3.5)] and will be extended 
rigorously in a forthcoming paper [TU] . Conjecture [1] has been proved in several cases: 
(a) for m = n < 9, using the exact values from [2]; (b) for sufficiently large m,n when 
st = o((mn) 1//5 ), using the asymptotics derived in [21]; (c) for several thousand values 
of (m, s; n, t) for m,n < 30. It has also been established to a high degree of confidence 
for many larger sizes using a simulation similar to the one described in [8], which was a 
variation on a method of Chen, Diaconis, Holmes and Liu [TT] . 

An interesting point of contrast between the integer and 0-1 cases is suggested by the 
above. Let E (m, s; n, t) be the quantity corresponding to E(m, s; n, t) in Remark [1] when 
the number of binary tables is decomposed in the same manner (see [7]). It was proved 
by Ordentlich and Roth [32] that E (m,s;n,t) < 1. A corollary of Corollary [T] would be 
that the opposite bound E(m, s;n,t) > 1 holds in the integer case. (In fact the smallest 
value we have found for s, t > is E(2, 3; 3, 2) = |||.) 

Some numerical examples comparing our estimates to other estimates and the correct 
values appear in Table [H 

Throughout the paper, the asymptotic notation 0(f(m,n)) refers to the passage of 
m and n to oo. Generally, the constant implied by the O may depend on a, b, and the 
sufficiently small and fixed constant e introduced later in the definition (13.11) of region 
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Table 1: Comparison of six estimates for M(m, s; n, t). 



m, s, n, t 


G(m, s; n, t) 


OH (B2.22)] 


m 






[2H Cor. 5.1] 


flH> 


Conjecture [T] 


Exact 


3,100,3,100 


7 

1.019x10 


1.192x10 


7 

1.262x10 






1.022x 10 


■1 /"'Of"! 1 A 1 

1.680x10 


(1.316±0.217)xl0 


13268976 


3,98,49,6 


7.594xl0 67 


3.716xl0 68 


1.278xl0 68 






9.630x10 


1.252x10 


(1.017±0.020)x 10 


l.OHOOx 10 


3 ^ 99; 9 j 33 


91 

2.116x10 


91 

2.788x10 


91 

2.864x10 






2.506x 10 


3.488x10 


(2.844±0.236)xl0 


2.79207x 10 


10,20,10,20 


7.434x10 


9.021 xlO 58 


1.511xl0 59 






1.059xl0 59 


1.226xl0 59 


(1.119±0.056)xl0 59 


1.09747xl0 59 


18,13,18,13 


5.157xl0 127 


6.962 xlO 127 


5.109xl0 130 






7.850 xlO 127 


8.502 xlO 127 


(8.065±0.224)xl0 127 


7.94500xl0 127 


30,3,30,3 


1.404xl0 92 


7.496 xlO 92 


1.086xl0 138 






2.223xl0 92 


2.315xl0 92 


(2.242±0.037)xl0 92 


2.22931 xlO 92 



TZ. We also use a modified notation 0(f(m,n)). A function g(m,n) belongs to this class 
provided that for some constant c, which is independent of a, b, and e, we have 

g(m,n) = 0{f{m,n)n C£ ). 

Under the hypothesis (11.11) of Theorem [H we have logm ~ logn, n = O(mlogm), 
m = O(nlogn), and A -1 = O(logn). Consequently, in the newly introduced notation, 
n = 0(m), m = 0(n), and A -1 = 0(1). In general, if ^,02,03,04 are constants, then 



2 An integral for M(m, s; n, t) 

We express M(m, s; n, t) as an integral in (m+n)-dimensional complex space then estimate 
its value by the saddle-point method. 

It is clear that M = M(m, s; n, t) is the coefficient of x{ ■ ■ ■ x s m yi • ■ • y n in 

m n 

n IK 1 -»**)"■ 

j=l k=l 

Applying Cauchy's Theorem we have 



/<->_ -\m+n r r s+1 s+1 t+1 t+ 

(2ttz) J J x x ■■■x m y 1 ■■■y n 



M = , n _^m+n f - f s+ l h s+1 t+l 1+1 dx l--- dx m ^1 ' ' ' dy n , (2.1) 
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where each contour circles the origin once in the anticlockwise direction. 



It will suffice to take the contours to be circles; specifically, we will put Xa = re t6j and 



y k = re k for each j, k, where 



A 



1 + A 
This gives 

M = (^^(^ A (l + A) 1+A ) m "/(m,n), (2.2) 

where 

n ^ r r n i . fc (i-A(e^ ) -i))- 1 ^^ 



where = (9 1 , . . . , 9 m ) and (p = . . . , n ). Let F(0, </>) denote the integrand in equa- 
tion f!2T3j) . 



3 Evaluating the integral 

This section follows closely the corresponding section in [7J. However, all proofs are 
intended to be complete, and not require the reader to refer to the latter work, except for 
the following lemma which we restate here without proof. 

Lemma 1. Let e, e" , s", e, A be constants such that < e < e" < e" , e > 0, and 
< A < 1 . The following is true if e" and e are sufficiently small. 

A A At f 

Let A = A(N) be a real-valued function such that A(N) = Q(N~ £ ). Let a,j = <L(iV) ; 
Bj = Bj(N), C jk = C jk (N), Ej = Ej(N), F jk = F jk (N) and jj = Jj(N) be complex- 
valued functions (1 < j,k < N) such that Bj, Cj k , Ej, Fj k = 0(N £ ), hj = 0(N 1 ^ 2+£ ), and 
Jj = 0{N~ l l 2+£ ), uniformly over 1 < j, k < N. Suppose that 

, N N N N 

f(z) = exp I -AN z j + Y + N Y Bj z j + Y CjkZjZk 
j=i i=i j=i j,k=i 

N N N s 

+ N E 44 + Y P ^ 2 A + J2 Jj z 3 + s ( z n 

3=1 j,k=l j=l ' 

is integrable for z = (z x , z 2 , ■ ■ ■ , z N ) G U N and S(N) = max z&u \5(z)\ = o{l), where 

U N = {z\ \ Zj \ < N~ 1/2+£ for l<j<N}, 

where i = e(N) satisfies e" < 2e < e" . Then, provided the 0() term in the following 
converges to zero, 

/ /(*)&=( ) exp(e i + e 2 + 0((N- A + 6(N))Z)), 
Ju N \AN J 
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where 



©i = — - — <ij h — ^ — o $ "i — ^ — ^ — ^ — o BjC jk 

2AN f-f 3 AA 2 N 2 J 16A 3 iV J A/" J 3 

3=1 3=1 3=1 3,k=l 

1 N 3 N 1 N 

+ w¥¥ CifeC ^ + il 2 ^ ^ Ej + ii^F ^ Fjk 

3,k,£=l 3=1 3,k=l 

1 N 3 ^ 45 ^ 

6A 3 N 3 ^! 3 2A 3 N 2 ^ 3 3 IQA 4 N 2 ^ 3 3 

]= 1 3=1 J= i 

1 N 3 N 1 N 

3,k=l 3=1 J,k=l 

1 N * ■ 3 N 

+ : " 54 , r 4 + 2a k )C jk C je + ~ 4 3 2 > (2d • + a k )BjC jk 

3,k,t=l 3,k=l 

/ i N N 

Z = exp — - V Im^) 2 + V Im(4-) 2 

3 N 1 N \ 

8A Ar ~ 16A d iV d -tTli / 



j,k=i i,M=i 

We use the notation 7£ c for the complement of a region 7?.. To evaluate the integral 
I(m, n) defined in (I2.3p . we proceed as follows: 

I(m,n) = [ F+ ! F= [ F + 0(1) [ \F\, K'DTZ. 
Jn Jn c Jn' Jn c 

The region TZ C [—ir, +n] rn+n is defined below in (I3.ip . The larger containing region TZ' 
is defined in (13. 7p . The asymptotic value of the integral J",/ F is obtained by substitut- 
ing TZ for the variable TZ" in equation (13.21) . and applying the main result of this section, 
Theorem [2j In the next section we show that for a certain constant c 5 > 

F\ = (3( e _c 5 min(m2e ' n2e )/ lo g ri ) / p 
tc Jn' 

(Again, the quantity e is a small positive constant arising in the definition (13.11) of the 

region TZ.) This is the complete summary of how we shall evaluate I(m,n), and now we 

may proceed to the technical details. 

For any region TZ we set 

I n "(m,n) = / F. 

Jn" 

In the integral (12.31) . it is convenient sometimes to think of 9j, <p k as points on the unit 
circle. We wish to define "averages" of the angles 8j, (fi k . To do this cleanly we make 
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the following definitions, as in [8]. Let C be the ring of real numbers modulo 2tt, which 
we can interpret as points on a circle in the usual way. Let z be the canonical mapping 
from C to the real interval (— it, it}. An open half circle is C t = (t — n/2, t + vr/2) C C for 
some t. Now define 

C = {x = (x 1 , . . . , x N ) G C N \x 1 , . . . ,x N G C t for some t G K}. 

If x = . . . , Xjv) G C then define 

AT 



More generally, if a; G then define a; = t + (x x — £, . . . , x N — t). The function a; i— > x 



is well-defined and continuous for x G C N . 



For some sufficiently small e > 0, let 7?. denote the set of vector pairs 0,4> G C xC 
such that 

|0 + 0| < (l + A) _1 (mn)- 1/2+2e 
region ft: \6j\ < (1 + A) -1 n~ 1/2+£ , 1 < j < m (3.1) 

|0 fc | < (1 + A)" 1 m~ 1/2+£ , 1 < k < n, 

where 9j = 9j — and <ft k = <p k — (j>. In this definition, values are considered in C . 

Let = (9 1 , . . . , m _i), <j> = . . . , n _i), and define T x to be the transformation 
Ti(9, 4>, fi, 5) = (0, 4>) given by 

H = + S = 0- c(>, 

together with 9j — 9j — 9 (1 < j < m — 1) and fc = fe — (1 < k < n — 1). We also 
define the 1-many transformation T* by 

Tl(0,&ri = \jT 1 (0,&^5). 

5 

After applying the transformation T\ to In{ m , n )) the new integrand is easily seen to be 
independent of 5, so we can multiply by the range of 5 and remove it as an independent 
variable. Therefore, we can continue with an (m+n—1) -dimensional integral over the 
region S defined by 7Z = 7i*(«S). More generally, if S" C [— |vr, \v;] m+n ~ 2 x [— 2tt, 2tt] and 
TZ" = T*(S'), we have 

I^ii(m,n) = 2Timn / G(0, (f), (i) d0d(f)dfi, (3-2) 
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where G(6, 0, ji) = F(Ti(0, 0, /i, 0)). The factor 2irmn combines the range of 5, which is 
47r, and the Jacobian of 7\, which is mn/2. 

Note that 5 is defined by virtually the same inequalities as define TZ. The first in- 
equality is now < (1 + A) 1 (mn) _1 ' 2+2E and the bounds on 

m— 1 re— 1 

m = 0j alld = ~ 

i=l fc=l 

still apply even though these are no longer variables of integration. 

Our main result in this section is the following. 

Theorem 2. Define A = |A(1 + A). Under the conditions of Theorem^ there is a region 
S' D S such that 

/ s ,G(9,0,rt^=M-/ 2 (_) (_) (_) 

Proof. For re real and |Ax| small, 

(1 - A(e ix ' - l))" 1 = exp(Xtx - Ax 2 - tA 3 x 3 + A 4 x 4 + 0((1 + A) 5 |x| 5 )) 

with 

A = |A(1 + A), A 3 = ±A(1 + A)(l + 2A), A 4 = iA(l + A)(l + 6A + 6A 2 ). 

We are about to state an estimate for G, and some further estimates fl 3 . 3 H in a moment, all 
of which hold uniformly in the region S. Let us alert the reader that we shall be defining 
a larger region, see (13.81) below, which contains S and in which all of these estimates 
continue to hold true. Uniformly in S, where all (1 + A)|/i + 9j + <j) k \ are small, 

G = exp j-A J^j > + h + <^) 2 - iA * J2j ^ + + ^) 3 

+ A 4 X) if > + h + + ((! + A ) 5 E j)fe ^ + + 4| 5 ) }■ 

Here and below, the undelimited summation over j, k runs over 1 < j < m, 1 < k < n, 
and we continue to use the abbreviations 9 m = — ^2^=1 @j an d — ~ J2k=i 4>k- 

We now proceed to a second change of variables, (0, <fi, /i) = T 2 (cr, r, /x) given by 

e j = (Tj + c/i 1; fc = r fc + dz/ x , 

where, for 1 < h < 4, /i h and v h denote the power sums Y^=i a j an d J2k=i r fc> respec- 
tively. The scalars c and d are chosen to eliminate the second-degree cross-terms cr„- aj and 
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r ki r k2 , and thus diagonalize the quadratic in a = (a ± , . . . , <7 m _i) and r = (r 1; . . . , T n _i). 
Suitable choices for c, d are 

1 ; , d 1 



1/2 ' 1/2 ' 

m + m ' n + n ' 



and we find the following: 

} + 6j + 0fc) 2 = rnn/j 2 + n/i 2 + 77M/ 2 

(/i + ^- + fc ) 3 = 3/i(n/i 2 + 77M/ 2 ) + ^(^3 + 3c/i 2 /ii) + ^(^3 + ?>dv 2 V\) 

+ 0((1 + A)-V 1/2 ) (3.3) 
(/i + 7 - + fc ) 4 = 6/i 2 z/ 2 + ^4 + mv A + C((l + A)~ 4 n~ 1/2 ) 
^ \ f , + § j + 4> k \ 5 ^d((l + X)- 5 n- 1 ^. 

The Jacobian of the matrix T 2 is (mn) -1 '' 2 , and so 

/ G = (mn)- 1/2 [ E x , (3.4) 

where E 1 = exp(L 2 + 0(n -1 ^ 2 )), and 



L 2 = —Arifi 2 — Amu 2 — Amnfi 2 — iA 3 nfi 3 — iA 3 mu 3 — 3iA 3 nfifi 2 — 3iA z m\iv 2 
SiA^cn^fix — ?>iA 3 dmv 2 v l + A 4 nfi 4 + A 4 mv 4 + 6A 4 fi 2 u 2 . 



(3.5) 



Define the regions Q, A4, S', and 1Z' by 

Q = { < 2(1 + A) _1 n _1/2+£ , j = 1, . . . , m-l } 

n{|r fe | < 2(l + A)~ 1 m" 1/2+£ , fc = l,...,n-l}, 
n{H <2(l + A)- 1 (mn)- 1/2+2e } 
M = { |^| < (1 + \y l m l/2 n~ 1,2+2£ } n { h| < (1 + A)-V /2 m- 1/2+2e }, 

S' = T 2 (QnM), (3.6) 

and 

K' = T?(S'). (3.7) 

Summing for 1 < j < m — 1 the equation 9j = <jj + c/i l5 and inserting the value of c, we 
find 

m— 1 

i=i 
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In the region S we have | Yl^=i @j\ — 0- + ^) ln +e , and so in T 2 1 (iS) we have 

IjtfjJ < (1 + A) m ' n ' . 
Using this, and the dual inequality for 1^1, the reader can check that 

T 2 \S) CQnM. 

It will be convenient if we can apply the bounds (13.31) throughout the expanded region 
S' = T 2 (<2n A4), rather than only in S. To see that this is valid, note that the calculations 
leading to these bounds hold equally well if the coefficient of e in an exponent of m or n is 
increased; or, if an assumption such as < (1 + A) -1 {mn) l ^ 2+2£ is made more permissive 
by a multiplicative constant: \fi\ < 2(1 + A) 1 (mii)'' 2+2£ . Therefore, it suffices to note 
that S' lies inside the region 

\H\ <2(l + \)- l (mny l/2+2e 
\0j\ < 3(1 + A) -1 n- 1/2+£ , 1 < j < m - 1 

\L\ < (I + ^ n- 1/2+2£ , (3.8) 
|4| < 3(1 + A)" 1 m~ 1/2+£ , 1 < k < n - 1 

Kl^i + A)- 1 ™- 1 / 2 ^. 

Define E 2 = exp(L 2 ). We have shown that in the region Q n M. the integrand E 1 
satisfies E l = E 2 (l + 0{n~ l/2 )). We can approximate our integral as follows: 

/ E 1= [ E 2 + d(n- 1 ' 2 ) [ \E 2 \ 
JQnM JQnM JQnM 

E 2 + 0{n~ l/2 ) I \E 2l 
QnM J Q 

(3.9) 



/ E 2 + 0(1) [ \E 2 \+0(n- 1/2 ) [ \E 2 \. 
JQ J QnM c JQ 



' Q JQnM JQ 

It suffices to estimate each of the three integrals in the last line of (13. 9p . 

We first compute the integral of E 2 over Q. We proceed in three stages, starting with 
integration with respect to fi. For the latter, we can use the formula 

/7T 1/2 / Q \ 

exp(-Amn/, 2 -^)^= {£-) ^— + 0{ n^ t 

— (mn) ' 

provided f3 = o{A{mn) l ^ 2+2e ). In our case, {3 = 3A 3 (n/j 2 + mu 2 ), which is small enough 
because m = 0(nlogn) and n = O(mlogm). Integration over \i contributes 

/ \ l ' 2 exp r -9At(nfi 2 + mu 2 ) 2 + ,-u\ (3 1Q) 

\AmnJ \ AAmn I 
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The second step is to integrate with respect to cr the integrand 

exp ( - An/j, 2 - j^lA - iA 3 nfi 3 - 3iA 3 cri[i x ii 2 

T x (3-H) 



This is accomplished by an appeal to Lemma [TJ We must take the N of that lemma 
equal to m — 1, the number of variables. This dictates that the limits of integration be 
±(m-l) 1/2+£ ~,but our limits, based on the definition of Q, are ±2(l + A)~V 1/2+£ . Thus, 
we make a change of scale, 



*3 4 



(l + A)(n/m) 1/2 

before integrating. This change of scale will introduce a multiplicative factor (1+A) ~^ m ~ 1 \m/n) 
into our evaluation of the integral. In the terminology of that lemma, we have N = m — 1, 
S(N) = 0(n~ 1 ), e = §£, e" = |e, e" = 3s, e = 6e, and e(N) = £ + o(l) is defined by 
2m~ 1 / 2+e — (m — l) _1 ^ 2+e . Furthermore, taking account of the scale, we find 

A m A 12AA 4 -9Alm 

A ~ 772 7' a i ~ „ „ /„ 772 ^2; 



;i + A) 2 m-l' 3 ~ 2.4(1 + A) 2 



n 



- iA 3 m 1//2 m ~ 3iA 3 m 3 ^ 2 

j = "(1+A) 3 n l/2 m-r ° 33 ' = "(TTAf 

~ A 4 m m ~ 9Aj m f 

&n = 7 1 F AA > = J , Ja = U. 

3 (1 + A) 4 nm-1' w 4,4(1 + A) 4 n' J 
We can take A = |, and calculate that 

Ej + > F../ = — — § + 0{n 



4i 2 A^ 



i =1 j,j'=i 



E 4 + E 44/ + — ^ — 3 E 4'4" = + 

1 V- 1 m(3A A 9Aj\ ~ _ x 

— — > cia H ^ — r > a,- = — — -y)p 2 + 0(n ) (3.12) 

2AN ~[ 3 AA N ~y 3 n\ A AA 2 J 

After checking that @ 2 = 0{n~ l ) = o(m~^ 5 Zi), we conclude that integration with 
respect to cr contributes a r-free factor 



, , , exp —I 1 _ + 0(m" 4/ %) . (3.13) 
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By the conditions of Theorem [TJ Z l < n 1 ^ 5 , so m *' 5 Z 1 = o(l) as required by Lemma [TJ 

Finally, we need to integrate over r. Collecting the remaining terms from (13.51) . and 

the terms involving r from (13.101) and (13.121) . we have an integrand equal to 

/ , /3A 4 m 9Alm\ 9 Aim 2 

e W ^-Amu 2 + - ^-)v 2 - _„ 2 

+ y4 4 mz/ 4 — iA 3 mu 3 — d>iA^dmv 2 Vi + (^(m -4 ^ 5 ^!^ 

Again we use LemmadJ This time, the factor due to scaling is (l+A)^ n ~ 1 - > (n/m) (ri ~ 1 ' ) / 2 , 
and N = n - 1, 5(JV) = C^m -475 ^), e' = §e, e" = §e, e"' = 3e, e = 4e, and 
e(iV) = e + o(l), as defined by 2n~ 1/2+£ = (n - l)~ 1/2+e *. The substitution table this 
time reads 

A n A 



^ ~~ T. 772 7> ^fc 



1 + Afn-l' 4A\1 + X) 



2 ' 



r ^3 n /2 n a 3iA 3 n 3/2 

fe " (l + A) 3 V/ 2 n-l' (1 + A) 3 dm m 3/2' 

fe (l + Afmn-l' fefc '~ 4A(1 + A) 4 m' " k ^ 

This time we take A 2 > 4/5 so that rT^ 2 = o(8(N)); calculations similar to the previous 
case lead to 



fc=1 fc,fc=i 
jv g jv 1 ^ 

16A 3 Nt^ 8A N V 16A 3 iV 3 ftT 

3A 2 n ~ _! 
+ On 



8A 3 m 

1 v^- 1 v^.2 3A 4 ~, _ u 

A AT ^— ' * </|2»r2^ * S/l 3 9 4 2 v ' 



2AN^ AA z N Zj ^ 8A* 2A Z 

= Z 2 = exp(^+d(n- 1 ))=0(l)exp 
\8A m J 



(1 + 2A) 2 n 
24Am 



Again 2 = 0(n 1 ), implying that 2 = o[m 4 ' 5 ZiZ 2 ). Including the contributions 
from (13.101) and (I3.13p . we obtain 



E 2 = (- 



7T \V 2 / 7T \(™~l)/2/ 7T xC™- 1 )/ 2 



q V ylmn / V An ) 



, 3A 4 (m n\/3A 4 15A 3 \ _ 4 / 5 

x exp o H 9 + 1 — n o + U(m ' Z X Z 2 

F V S4 3 24 2 \n m) \AA 164 3 / V 1 2 ' 



(3.14) 
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Since 

Z X Z 2 = 0(1) exp ( ^ ( — + — 
• 1 8A 3 V n m 



3A^ (vti n 



and since 

3Aj _ (1 + 2A) 2 
8^ ~ 2AA ' 

it follows from the main hypothesis (11. ip of Theorem [1] that Z X Z 2 = 0(n 6a ^ 5 ). The 
condition a + b < 1/2 implies — 4/5 + 6a/5 < —6 — 1/5; hence, substituting the values of 
A, A 3 , A 4 , we conclude that 



/ 7T \!/ 2 / 7T \('«-l)/2/ 7T \(l-l)/2 

2 Vv4mri/ VAn/ Vv4m/ 

3.15 

x exp — — + — + 0(n b 1/5 

We next infer a estimate of \E 2 \. The calculation that lead to (I3.14p remains 
valid if we set A 3 to zero, which is the same as replacing L 2 by its real part. Since 
\E 2 \ — exp (Re(L 2 )), this gives 

f .„ . /(1 + 2A) 2 / 5n 5m\ .A f „ 

= 0(n°) / £ 2 (3.16) 
J Q 

under the assumptions of Theorem [TJ The third term of (13.91) can now be identified: 

0(n- 1/2 ) [ \E 2 \ = 0(n- 1/2+a ) [ E 2 = 0(n h ) [ E 2 . (3.17) 



Finally, we consider the second term of (13. 9p . namely 

\E 2 \, 

QnM c 

which we will bound as a fraction of J Q \E 2 \ using a statistical technique. The following 
is a well-known result of Hoeffding 



Lemma 2. Let X X ,X 2 , . . . ,X N be independent random variables such that EX i = and 
— M for all i. Then, for any t > 0, 

Prob ( V X > t] < exp ( t —s 

\^ ~ J ~ V 2NM 

\ l= i / \ 
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Now consider \E 2 \ = exp(Re(L 2 )) . Write M = M 1 H M 2 , where M 1 = { \iix\ < 
m 1//2 n _1//2+2£ } and Ai 2 = { |^i| < n 1//2 m~ 1 ^ 2+2£ }. For fixed values of /i and <x, Re(L 2 ) 
separates over r 1; r 2 , . . . , r n-1 and therefore, apart from normalization, it is the joint den- 
sity of independent random variables X 1 ,X 2 , . . . , X n _± which satisfy EX k = (by sym- 
metry) and |Jf fe | < 2(1 + A) _1 m _1 ^ 2+e (by the definition of Q). By Lemma [21 the fraction 
of the integral over r (for fixed li, <t) that has v x > (1 + \)^ 1 n l ^ 2 mr 1 ^ 2+2e is at most 
exp(— m 2e /8). By symmetry, the same bound holds for v 1 < —(1 + X)^ 1 n 1 ^ 2 m^ 1 ^ 2+2£ . 
Since these bounds are independent of ll and <r, we have 



/ \E 2 \ < 2exp(-m 2e /8) / \E 2 
jQnM c 2 Jq 



By the same argument, 



\E 2 \ < 2exp(-ra E /8) / \E 2 \. 
QnMl JQ 



Therefore we have in total that 



/ |E 2 | < 2(exp(-m 2 78) +exp(-n 2 78)) / \E 2 \ 
Jqhm c Jq 

< 0(n- b ) [ E 2 , (3.18) 
Jq 

using again (13151) . Applying (j33j) with fl3~T5l) . f[3TTTj) and (13TT8D . we find that j QnM E 1 
is given by (13.151) with the error term replaced by 0{n~ h ). Multiplying by the Jacobian 
of the transformation T 2 , we find that Theorem [2] is proved for S' given by (13.61) . □ 



4 Concentration of the integral 

Recall that E(0, 4>) is the integrand in equation (12.31) defining J(m, n), and that 1Z is the 
region given by (13.11) . In the previous section we estimated the integral of E(0, cf)) over a 
particular superset TZ' ~D 71. In this section we show that the integral of F(0, <fi) outside 
1Z is negligible in comparison if A is polynomially bounded. Larger values of A will be 
handled in the following section. 

Theorem 3. Suppose that m, n — > oo in such a way that (11. lft holds and A = n ^. 
Define I by 

v ; \AmnJ \An) \AmJ \ 24 A \n m) J K ' 

Then, for sufficiently small e > 0, 

/ \F\ = O(e- nS )I . 
Jn c 
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We begin with two technical lemmas whose proofs are omitted. 
Lemma 3. The absolute value of the integrand of I(m,n) is 

1^(0,0)1 =n/(*i+&), 

where 

f(z) = (1 + 4A(1 -cos^))" 1/2 . 
Moreover, for all real z with \z\ < -^(1 + A) _1 ; 

0< f{z) <exp(-Az 2 + {±A + A 2 )z 4 ). □ 

Lemma 4. Define N = [6000(1 + A)] , 5 = 2tt/N and g(x) = -Ax 2 + (\A + 27A 2 )x 4 . 
Then, uniformly for X > 0, K > 1, 

/30S 
exp(Kg(x)) dx < y/n/(AK) exp(0(Ar _1 + (AK)~ 1 )). □ 
-305 

Proof of Theorem^ Let iV and 5 be as given in Lemma 4. Define the region A to be the 
set of those (0, (f>) such that 

cos(9j + 4> k ) < cos 5 

for at least | min(mn e , m £ n) pairs (j, k). Define x ,x x , . . . , x N _i by x e = £5. 

If X C (— 77, 7r], we denote by N (X) the number of values of j such that 9j G X, 
and similarly define N^(X). Define region TZi(£) to be the set of those (0,0) such that 
N g ([x e - 45, x e + 45]) >m-m £ and N^-Xf, - 45, -x e + 45]) >n-n £ . 

Let U = U^tj 1 ^i (^) • The P ro °f °f t ne theorem consists in proving these three rela- 
tions: 

AUU = [-7T,7r] m+n (4.2) 
/ |F| = O(e- n )lo (4.3) 
/ |F| = 0(e-»")J a . (4.4) 

To begin the proof of (14.21) . we show that any point (6, <fi) for which A^ ([^ — 5, x £ +5] > 
m £ belongs to A U U. Indeed, if such a point does not belong to A, then it must have 
N ( f ) ([—x e — 25,—x £ + 25}) > |n. This in turn forces N g ([x e — 35, x £ + 35]) >m — m £ , which 
forces A^Q— x e — 45, — + 45]) >n — n . In particular, (0, <fi) e 7£i(£). 

To complete the proof of (14.21) . we show that (0, 0) belongs to A if A^Qa^ — 5, x e +5]) < 
m £ for all I. Let a be a minimum-length interval [a^,ay] such that N 9 (a) > m £ . Then 
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Nq(cl) < 2m £ and the complementary interval a has N (a) > m — 2m £ . If b = [x e > +2 , x e _ 2 ] 
(a subinterval of a), then N e (b) > m — 4m e . Thus, there are at least m £ disjoint pairs 
(j, j') with 9a G a and 9-> G b. (By disjoint we mean there are 2m £ distinct indices j 
involved.) Because of the 25 spaces between arcs a and b, for each k and each pair (j, j') 
at least one of cos(0j + <f) k ) or cos(0y + 4>k) is bounded above by cos 5. This implies that 
(0, <f>) G .A, as claimed, and completes the proof of ( 14. 211 . 

We turn next to ( 14.31) . Since A5 2 = @(A(1 + A) -1 ), Lemma 3 implies 

\F(6, <p)\ < exp(— CxA(l + A) -1 min(m £ n, mn £ )) 

for all (0, <fi) G A and some c x > 0. The volume of A is no more than (27r) m+n , and 
A(l + A)" 1 > (hgn)' 1 by (HID, so 

/ \F\ < (27r) m+n exp(— c x (log n) _1 min(m £ n, mn £ )) . 

From (11.11) . which implies that A = O(logn), and the assumption that A = n ^ 1 , we have 
that Jq = exp(0(mlogn + nlogm)). Relation 04.3(1 follows. 

Finally we come to (jO]l . For (0» G ft^) we define 5 = S o (0), Si = S x (0) and 
5*2 = S 2 (0) to be the set of indices j such that \9j — x e \ is: less than or equal to 45, in the 
interval (45,55], and larger than 55, respectively. The index sets T = T (</>), 7\ = Xi(</>) 
and T 2 = T 2 {4>) are defined similarly. Define lZ 1 (£;m 2 ,n 2 ) to be that subregion of TZx(£) 
for which |S 2 | = m 2 and |T 2 | = n 2 . Define U(m 2 ,n 2 ) by 

U(m 2 ,n 2 ) = \^jn 1 (£;m 2 ,n 2 ). 

e=o 

We note that m 2 and n 2 vary over the ranges [0, m £ ] and [0,n e ], respectively. Define 
U Q = W(0, 0) and U*=U\ U . 

Suppose (0, <p) G 72-i (0; m 2 , n 2 ). If \9j\ and |0 fc | are both less than or equal to 55, then 
Lemma 3 is applicable to f(9j + <p k ). If one of the two is less than or equal to 45 and 
the other exceeds 55, then cos(6 ) J - + <p k ) < cos 5. Thus, there exists c 2 > such that for 
(0,0) G^ 1 (0;m 2 ,n 2 ) 

' exp(-A(B j + fe ) 2 + (i^ + A 2 ){9 3 + fc ) 4 ), 

if (j^je^u^x^uTO, 

(4.5) 

exp(-2c 2 A(l+A)- 1 ), if (J,k) G (S x T 2 ) U (S 2 x T ), 

1, otherwise. 



m+<Pk)< { 



Since (0, 0) G U, we have |S | >m — m £ and |T | > ri — n £ ; thus, the size of (S x T 2 ) U 
(5 2 x T ) exceeds \{mn 2 + nm 2 ). Integrating the upper bound on |F| implied by (14. 51) . 
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we find 



In 



F\ < (27r) m ^( m )( n ) exp(-c 2 A(l+A)- 1 (mn 2 + m 2 n)) l' 2 {m 2 ,n 2 ) (4.6) 

) \ m 2/ \ n 2/ 



'7J 1 (0;m 2 : n 2) 

with 

/5S r55 / 11 II \ 

■■■ exp (0, + <t>k? + ih A + A 2 ) {O j + fc ) 4 dd"d<P", 

-58 J-58 V ' 

where the double-primes denote restriction to j G S U Si and k G T U T x . The factor 
(27r) m2+n2 comes from integrating over 9j for j G 5* 2 and 4> k for /c G T 2 , while the binomial 
coefficients account for the choices of S 2 and T 2 . 

Set m = m — m 2 , n — n — n 2 . As implied by the notation, the value of I 2 (m 2 , n 2 ) is 
independent of which specific variables constitute the sets S U Si and T U T x . 

Summing on m 2 + n 2 > 1 yields an upper bound for the integral of |F| over the 
region Um 2 +ri 2 >i m 2? ^2)- Notice, however, that the transformation 9j \— > — a^, 
0fc l—> 0fe + ^ is a volume-preserving bijection of 7Zi(£;m 2 ,n 2 ) with 7£i(0; m 2 , n 2 ) which 
leaves F(0, <f>) invariant. Thus, introducing an additional factor of iV to account for all 
values of £, we have 

, if! ^^(:)(:) (2 ^ + " 2 

x exp(— c 2 A(l + A) l {mn 2 + nm 2 )} l' 2 (m 2) n 2 ). (4.7) 

We now analyze / 2 (m 2 ,n 2 ) more closely. Define = (m')~ £ 6j, 9j = 9j - 
for j G S y S u (j) = (n'y 1 E'Vfc; <Pk = <Pk~ 4>' for k G T U 7\, // = <// + 6' and 
v' = 4> — 6 . In terms of //, z/', fc we have 

E ft + 

= m'n'fi' 2 + n ^ . Q) + m' 

and 



E + ^ 27mV/i' 4 + 27n' ]T . 6> 4 + 27m' ^ 4 . 



-f- zin y 

The latter follows from the inequality (x + y + z) A < 27(x 4 + y 4 + z 4 ) valid for all x, y, z. 
It follows that 

•• y e^\m'n'g(fj!) + n '^2 j 9(Pj) + m ' '^2 k 9(4>k)) dO"d<j>". (4.8) 

For the moment we are thinking of //, u, 9j, and fe (1 < j < m , 1 < k < n) as functions 
of the variables of integration. Arbitrarily choose one element G S USi and one element 
K G T U Ti. Define 5 3 = (S U S^) - {j,} and T 3 = (T U 7\) - {A;,}. We claim that the 
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inequality (14.81) remains valid when ^\ = Y2jeS us 1 1S replaced by Yl'j^j = ^2jeS > anc ^ 
similarly for k*. To verify this assertion, note that since \9j\ < 55 for all j G 5 U it 
must be the case that \6'\ < 55, too, and so \9jJ < 105. The claim now follows because, 
as is readily checked, g(x) < for \x\ < 105. 

As in the previous section, the next step is to make a change of variables to //, v ', 9j 
(j G S s ), and fc (k G T 3 ). The region corresponding under this change of variables to 
[—55, 55] m n is not exactly a product, but it is contained in the product [—105, 105} m n . 
(The argument given a few lines earlier to show \9jJ < 105 applies equally well to \9j\, 
j G S 3 .) Hence, 

/10<5 plOS . 
■ ■■ exjp(m'n'g(fi') + n V] g(9j) 
10<5 J-108 V ^i^Sz 

+ rri ^2 keTs 9(4>S) d ^jes 3 d 4>keT s drf . 

Here, the factor 0(5)m'n comes from the integration of v over [—105, 105] and the Jaco- 
bian m'n/2 of the transformation. Now the integral splits into a product of m + n — 1 
one-dimensional integrals. Each of these factors can be bounded by Lemma HI and we 
find that 

^K,n 2 )<0(5)mV( I £ y ) 1/2 (-£ 7 ) (n (-£-,) ^ 
x exp(o(l) + 0(1 + A- x )(m'/ri + n'/m)). 

By m ~ m, n ~ n and (II. ip . the exp(- • •) term above equals n 0<yl \ Also, 

\Am'J " UrJ 

and similarly for the other terms. Noting that 5N = 0(1), we find from (14 .7\\ that 

[ \F\< ^° (m2+n2) exp(-c 3 A(l + A)- 1 (mn 2 + nm 2 ))/ . 

* m 2 +n 2 >l 

Since A = ^((logn) -1 ), the sum is dominated by the terms with m 2 + n 2 = 1. We 
conclude that 

' \F\ = O(e-" e )/ . (4.9) 



It remains to consider U D 72. c . As was argued before in obtaining (14.61) we have in 
the m 2 = n 2 = case 

/ \F\< [ exp(-Aj2(0 j + <Pk) 2 + (M + A2 )J2^ + ^ 4 ) ded ^ 
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where now the sums are over all j, k. Observe that the transformation 9j h- > 9j — x^, <p k h- > 
fc + takes T^^O, 0) D 7£ c bijectively to 72. 1 (0;0, 0) fl 7Z C , since 1Z,1Z C are invariant 
under this mapping. Hence, the integral over U fl TZ C of |F| is no larger than N times 
the integral on the right side of the last displayed relation. 

The region of integration, [—56, 56} m+n fl 1Z C , is contained in the union U/S"^ ^hi 
where V h equals the product [—56, 56} m+n intersected with that part of 1Z C where the hth 
inequality of definition (13. ip fails. As was done in the case, we may bound L • • • by a 
product of m+n — 1 one-dimensional integrals. However, in the present situation, when we 
make the change of variables to (//, v , 8j, <f) k ), the transformed region of integration, albeit 
still contained in the product [—105, 10<5] m+n ~ 1 (we are omitting v from the discussion), 
has the additional property that one of the latter intervals, the one corresponding to the 
index h, is missing a neighborhood of size (1 + A) _1 (m?2) _1 ^ 2+2e , (1 + A) _1 n 2+e , or 
(1 + A) _1 m" 1/2+£ about 0. 

Throughout [—105, 105] we have g(x) < —Ax 2 /2. Because of the inequality 

I" e~ KA * l2 dx < — ^— ex P (-±A(l + A)- 1 ^), 

J(i+X)-'k-^ ~ (KA) 1/2 V 5 ; 

(K being one of mn, m, or n; and a being, respectively, 2e, e, or e), one of the m + n — 1 
factors is smaller than (n / Amn) 1 ^ 2 , (it / Am) 1 / 2 , or (it/An) 1 / 2 (respectively, depending 
on K) times exp(— |A(1 + A) _1 min(m 2e , n 2e )) . (One may assure that neither of the indices 
j*, chosen for omission in the change of variables coincides with the index connected 
to the hth inequality of the negation of definition (13.11) .) The factor exp(— | • • ■ ) is, of 
course, independent of h. Allowing iV values of i and m + n + 1 values of h, we have 

f , , , nm / 7T \V2 / 7T \(«-l)/2 / 7T \ (m-l)/2 

Ln' m - {m + n + 1)N0(S)n (AST) (AS) (a) 
x exp(— gA(l + A) -1 min(m 2e , n 2e )) 

< e I . 

This inequality and (14.91) together imply (14.41) . thus completing the proof of Theorem [2j 

□ 



5 Proof of Theorem H 

Once again we remind the reader that F(6, </>) denotes the integrand in equation (12.31) 
defining I(m,n). We continue to use J for the quantity defined in equation (14.11) . 
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Proof of Theorem^ Let a,b be given positive numbers satisfying a + b < |, and let 
(m, n, s, t) be a sequence of 4-tuples such that m, n — > oo in such a way that hypothe- 
sis (II. ip is satisfied. 

We first note that both (11.31) and Corollary [1] are equivalent to ( 11.2jl for the respective 
ranges of A indicated. Consequently, to prove Theorem [1] and Corollary [1] it suffices to 
prove (11. 3p for A = 0(n 5 ) and Corollary [1] for A > n 5 . 

We will divide the proof into two ranges of the parameter A. First we assume that 
A = 0(n 5 ), where Theorem [3] can be applied, then we use a different method for A > n 5 . 

First suppose that A = 0(n 5 ). As explained in Section [3], 

M(m, «; n, t) = — * (A~ A (1 + A) 1+A p ( jf # F + 0(1) ^ |F|) , (5.1) 

where 1Z is defined by (13.11) and 1Z D 1Z is defined by (13. 7p . By (13. 2p and Theorem [21 the 
first integral on the right side of (15.11) equals 27rexp(| + O(n~ b ))l , and by Theorem [31 
the second equals O(n _1 )/ . These values yield Theorem [TJ for this case in the form given 
by (Q. 

For the remainder of the proof, we assume that A > n 5 . By the first part of this proof, 
Corollary [T] holds for A = 0(n 5 ); hence, 

MK2nV,2mn 5 ) =2Kl)M] (l + Q(n*)). (5.2) 
M(m, n ; n, mn ) 

Let s = lcm(m, n)/m, t = lcm(m, n) / n, and let V be the convex polytope of m x n real, 
nonnegative matrices whose rows sum to s and whose columns sum to t . For every pair 
(s, t) such that ms = nt there is an integer q such that (s, t) = (qs Q , qt ), and M(m, s; n, t) 
equals the number of integer lattice points in the dilated polytope qV. The latter count, 
as a function of the positive integer q, is called the polytope enumerator for V, and is 
denoted L v (q), pQ. Thus, we have 

M(m, s; n, t) = L v (q). 

The polytope V, an example of a transportation polytope, has integral vertices [6], and so 
Lp(q) is a polynomial in q, the Ehrhart polynomial. The degree of L v is d, the dimension of 
the polytope V\ in our case, d = (m— l)(n— 1). By a theorem of Stanley (an algebraic proof 
first appeared as Proposition 4.5 in [37]; a more geometric proof is given as Theorem 2.1 
in [36]) the representation of L v (q) in a particular basis, 



n V / 
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has all its coefficients h , . . . , h d nonnegative. For q > d, we have the expansion 



and so 



d ) V Q—j) \dj ' TTocv(<7— 7 

w^-^rw?-;) 



where g(fc) = J2i=o ^d-i^ > 0. Note that "P, and therefore g(k) for each k, depend only 
on m and n. Also note that q/ Yljexil ~ j) ls decreasing as a function of q for q > d and 
any non-empty X C [0, d— 1]. 

Since g = Agcd(m, n), we conclude that there is a function a(m, n, A) such that 

M(m, An; n, Am) = Q^^^g(0)(l + a(m,n, X)/\) (5.3) 

a(m, n, A) > for q > d (5.4) 
a(m, n, A) is decreasing in A for fixed m, n and q > d. (5-5) 

For A > n 5 , since d = 0(n 2 ), 
Hence, by (15.31) 

M(m,2n 6 ;n,2mn 5 ) _ (m _i )(w -i) , pc. _ M 1 + |a(ro, w, 2n 5 )/n 5 
M[m,n ;n,mn ) 1 + a(m, n, n )/n 

Comparing this to (15.21) and noting from (I5.4[l5.5p that < a(m,n, 2n 5 ) < ct(m, n,n ), 
we conclude that a(m, n, n 5 ) = 0(n 5 ~ b ). This implies by (I5.5P that a(m, n, A) = 0(n 5 ~ b ) 
for A > n 5 . Using this information about a(m, n, A) with (I5.3P gives 



M(m, An; n, Am) / A 



(m-l)(n-l) 



u , - m , (l + 0(n- 6 )). 

M{m,n -,71,17171 ) \n J 

Combining this with Corollary [1] for A = n 5 shows that Corollary [1] holds for all A > n 5 . □ 
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